EIGENVALUES OF VT- SYMMETRIC OSCILLATORS WITH 
POLYNOMIAL POTENTIALS 



KWANG C. SHIN 

Abstract. We study the eigenvalue problem —u"{z) — [{iz) m + P m -i(iz)]u(z) — \u(z) 
with the boundary conditions that u(z) decays to zero as z tends to infinity along the rays 
argz = — | ± where P m _i(z) = a\z m ~ x + a 2 z m ~ 2 + ■ ■ • + a m _iz is a polynomial and 
integers m > 3. We provide an asymptotic expansion of the eigenvalues A n as n — ► +00, 
and prove that for each real polynomial P m -i, the eigenvalues are all real and positive, with 
only finitely many exceptions. 
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1. Introduction 

For integers m > 3 fixed, we are considering the "non-standard" non-self-adjoint eigenvalue 
problems 

d 2 



[l) Hu(z,X):-- 



u(z, A) = Xu(z, A), for some A G C, 



alz 2 

with the boundary condition that 

7T 27T 

(2) u(z, A) — > exponentially, as z — > 00 along the two rays arg(z) = ± 



2 m + 2 

where P m _i is a polynomial of degree at most m — 1 of the form 

(3) P m _i(z) = ai^ m_1 + a 2 z m-2 H h a m _i^, aj G C for 1 < j < m - 1. 

We let 

a := (ai, a 2 , . . . , a m _ x ) G C m ~ l 

be the coefficient vector of P m _i(^). We are mainly interested in the case when P m -\ is 
real, that is, when a G R m_1 . However, some interesting facts in this paper hold also for 
a G C m_1 . So except for Theorem H] below, we will use a G C m_1 . 

If a nonconstant function u satisfies with some A G C and the boundary condition 
0, then we call A an eigenvalue of H and u an eigenfunction of H associated with the 
eigenvalue A. Also, the geometric multiplicity of an eigenvalue A is the number of linearly 
independent eigenf unctions associated with the eigenvalue A. The operator H in (0) with 
potential V(z) = —{iz) m — P m _i{iz) is called VT -symmetric if V(—z) = V(z), z G C. Note 
that V(z) = —(iz) m — P m _i(zz) is a PT-symmetric potential if and only if a G lR m_1 . 
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Before we state our main theorems, we first introduce some known facts by Sibuya 
about the eigenvalues A of H. 

Theorem 1. The eigenvalues X n of H have the following properties. 

(I) The set of all eigenvalues is a discrete set in C. 
(II) The geometric multiplicity of every eigenvalue is one. 

(III) Infinitely many eigenvalues, accumulating at infinity, exist. 

(IV) The eigenvalues have the following asymptotic expansion 

2m 

fr(^ + -) /7c(n--)\ m+2 

(4) X n = V2 . "V V . V . . 2j [1 + o(l)l as n tends to infinity, n G N, 

V Sin (^) F ( 1+ m) / 

where the error term o(l) could be complex-valued. 

This paper is organized as follows. In Section [21 we will introduce work of Hille ^2] and 
Sibuya [TH], regarding properties of solutions of (0). We then improve on the asymptotics of 
a certain function in [IB]- In Sectional we introduce an entire function C(a, A) whose zeros 
are the eigenvalues of H, due to Sibuya [TH]. In Section H] we then provide asymptotics of 
C(a, A) as A — > oo in the complex plane, improving the asymptotics of C(a, A) in [T%| . In 
Section we will improve the asymptotic expansion (0J) of the eigenvalues. In particular, 
we will prove the following. Throughout this paper, we use that [^J is the largest integer 
that is less than or equal to x G R. 

Theorem 2. Let a G C m_1 be fixed. Then there exist e^(a) G C, 1 < £ < y + 1 such that 
the eigenvalues X n of H have the asymptotic expansion 

Lf+iJ 

4^ i— l - ( i— 

(5) A n = A ,n+ V e t (a)X m +o[\l n 

n^+oo ' * \ 



where 



2m 
m+2 



V = ( } U + . 4Z\ ) Witk K m= I i Vl+t m - Vt™) dt > 0. 



One can compute K m directly (or see equation (2.22) in [0] with the identity T(s)r(l — s) 
7rcsc(7rs)) and obtains 

v^r(i + a 



Km 2cos(^)r(f + 

V in I V 2 ml 

In the last section, we prove the following theorem, regarding monotonicity of |A n |. 
Theorem 3. For each a G C' m_1 there exists M > such that \X n \ < |A n+ i| if n > M. 
This is a consequence of (JSJ). 

Finally, when H is PT-symmetric (i.e., a G M m_1 ), u(z, A) is an eigenfunction associated 
with an eigenvalue A if and only if u(—z, A) is an eigenfunction associated with the eigenvalue 
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A. Thus, the eigenvalues either appear in complex conjugate pairs, or else are real. So 
Theorem El implies the following. 

Theorem 4. Suppose that a G M m_1 . Then the eigenvalues A of H are all real and positive, 
with only finitely many exceptions. 

For the rest of the Introduction, we will mention a brief history of problem ((TJ . 

In recent years, these PT-symmetric operators have gathered considerable attention, be- 
cause ample numerical and asymptotic studies suggest that many of such operators have 
real eigenvalues only even though they are not self-adjoint. In particular, the differential 
operators H with some polynomial potential V and with the boundary condition (J2J) have 
been considered by Bessis and Zinn- Justin (not in print), Bender and Boettcher [2] and many 
other physicists [SI El IHl UHl U^l USl EIZI [THJ - 

Around 1992 Bessis and Zinn- Justin (not in print) conjectured that when V(z) = iz 3 + 
(3z 2 , /3 G M, the eigenvalues are all real and positive, and in 1998, Bender and Boettcher 
|2j conjectured that when V(z) = {iz) m + (3z 2 , j3 G R, the eigenvalues are all real and 
positive. Many numerical, asymptotic and analytic studies support these conjectures (see, 
e. g., [3 il U d 13 113 EE3 EEH] and references therein and below). 

The first rigorous proof of reality and positivity of the eigenvalues of some non-self-adjoint 
H in (PJ) was given by Dorey, Dunning and Tateo [B] in 2001. They proved that the eigen- 
values of H with the potential V(z) = — {iz) 2m — a^z)" 1-1 + ^±^i ; m, a,£ E K, are all real 
if m > 1 and a < m + 1 + \2£ + 1\, and positive if m > 1 and a < m + 1 — \2£ + 1\. 

Then in 2002 the present author ^H] extended the polynomial potential results of Dorey, 
Dunning and Tateo to more general polynomial cases, by adapting the method in [Hj. Namely, 
when V(z) = —{iz) m — P m _i(zz), the eigenvalues are all real and positive, provided that for 
some 1 < j < y the coefficients of the real polynomial P m ~i satisfy (j — k)ak > for all 
1 < k < m — 1. 

However, there are some PT-symmetric polynomial potentials that produce non-real eigen- 
values. Delabaere and Pham 0, and Delabaere and Trinh studied the potential iz 3 + ^iz 
and showed that a pair of non-real eigenvalues develops for large negative 7. Moreover, 
Handy [HJj, and Handy, Khan, Wang and Tymczak pi] showed that the same potential ad- 
mits a pair of non-real eigenvalues for small negative values of 7 ~ —3.0. Also, Bender, Berry, 
Meisinger, Savage and Simsek [T| considered the problem with the potential V(z) = z 4 + iAz, 
A G M, under decaying boundary conditions at both ends of the real axis, and their numeri- 
cal study showed that more and more non-real eigenvalues develop as \A\ — * 00. So without 
any restrictions on the coefficients a^, Theorem 0] is the most general result one can expect 
about reality of eigenvalues. 
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Figure 1. The Stokes sectors for m = 3. The dashed rays represent aigz = ±|, 7r. 

Also, the method used to prove Theorem 0] in this paper is new. The method used in 
[HI E| is useful in proving reality of all eigenvalues, but I think that some critical arguments 
in proving reality of eigenvalues in [HJ cannot be applied to the cases when some non-real 
eigenvalues exist. The asymptotic expansion (0) itself is interesting, and also © implies 
Theorem El Note that is not enough to conclude Theorem 01 Finally, Theorem El and 
PT-symmetry of H explained right before Theorem 0] above imply the partial reality of the 
eigenvalues in Theorem 



2. Properties of the solutions 

In this section, we introduce work of Hille [12] and Sibuya [TH] about properties of the 
solutions of ((H). 

First, we scale equation (JTJ) because many facts that we need later are stated for the scaled 
equation. Let u be a solution of (JH) and let v(z, A) = u(—iz, A). Then v solves 

(6) -v"(z, A) + [z m + P m ^{z) + AM*, A) = 0, 

where m > 3 and P m ~\ is a polynomial (possibly, P m _i = 0) of the form Q. 

Since we scaled the argument of u, we must rotate the boundary conditions. We state 
them in a more general context by using the following definition. 

Definition. The Stokes sectors Sk of the equation © are 

S k = LeC: 

See Figure ^ It is known from Hille [T2"l §7.4] that every nonconstant solution of © 
either decays to zero or blows up exponentially, in each Stokes sector S k - That is, one has 
the following result. 



aw) 



2kn 
m + 2 



< 



TV 



m 



for k e Z. 
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Lemma 5 ([TJ §7.4]). 

(i) For each k G 7L, every solution v of (with no boundary conditions imposed) is 
asymptotic to 



(7) (const.)z 4 exp 



as z — > oo m every closed subsector of Sk- 
(ii) // a nonconstant solution v of (J0J) decays in Sk, it must blow up in Sk-i U 
However, when v blows up in Sk, v need not be decaying in Sk-i or in Sk+i- 

Lemma El (i) implies that if v decays along one ray in Sk, then it decays along all rays in 
Sk- Also, if v blows up along one ray in Sk, then it blows up along all rays in Sk- Thus, 
since the rotation z t— > iz maps the two rays in (J2J) onto the center rays of S_i and Si, 

the boundary conditions on u in ([TJ) mean that v decays in S-i U Si. 

Next we will introduce Sibuya's results, but first we define a sequence of complex numbers 
bj in terms of the and A, as follows. For A G C fixed, we expand 

(1 + aiz- 1 + a 2 z~ 2 + ■■■ + a m _ lZ l - m + \z~ m ) l/2 

oo /u 

= i + X) ( I) ( ai2_1 + G2Z ~ 2 + • " + a — ^ 1_m + A ^ m ) fe 

fe=l ^ ' 



oc 



(8) = i + ^M^l, for large |z|. 

Note that 6 2) &m-i do not depend on A, so we write bj(a) = bj(a,\) for j = 
1, 2, . . . , m—1. So the above expansion without the Az _m term still gives bj for 1 < j < m—l. 
We further define r m = — -j if m is odd, and r m = — ^ — &m +1 (a) if m is even. 

The following theorem is a special case of Theorems 6.1, 7.2, 19.1 and 20.1 of Sibuya JH] 
that is the main ingredient of the proofs of the main results in this paper. 

Theorem 6. Equation with a G C m_1 , admits a solution f(z,a,X) with the following 
properties. 

(i) f(z,a,X) is an entire function of z,a and A. 

d_ 

dz • 



[ii) f(z, a, A) and f'(z, a, A) = jrf{z, a, A) admit the following asymptotic expansions. 



Let e > 0. Then 



f(z, a, A) = + O^z- 1 ' 2 )) exp [-F{z, a, A)] , 

f(z, a, A) = - (1 + Oiz- 1 ' 2 )) exp a, A)] , 



as z tends to infinity in the sector | argz| < — £ > uniformly on each compact set 
of (a, X) -values . Here 

F(z,a,X) = -^—z^ +1 + V bAa)z^ m+2 ^\ 

(iii) Properties (i) and (ii) uniquely determine the solution f(z, a, X) of (0J). 

(iv) For each fixed a G C" 1 " 1 and 5 > 0, f and f also admit the asymptotic expansions, 

(9) /(O, a, A) =[1 + o(l)]A- 1 / 4 exp [L(a, A)] , 

(10) f (0, a, A) = - [1 + o(l)]A 1/4 exp [L(a, A)] , 
as A — ► oo in the sector | arg(A)| < n — 5, where 

/ ( + Pm-i{t) + A - ^ - J2j=i bjia)^) dt if m is odd, 
jo +0 ° (V* m + Pm-i(*) + A - if - bjia)^ - ^g^) ^ i/ m even. 

(v) TTie entire functions A i— > /(0, a, A) and A h- > /'(0, a, A) nai>e orders ~ + — . 



L{a,X) 



Proof. In Sibuya's book |18j . see Theorem 6.1 for a proof of (i) and (ii); Theorem 7.2 for a 
proof of (iii); and Theorem 19.1 for a proof of (iv). Moreover, (v) is a consequence of (iv) 
along with Theorem 20.1. Note that properties (i), (ii) and (iv) are summarized on pages 
112-113 of Sibuya PH- □ 



Using this theorem, Sibuya [T%1 Theorem 19.1] also showed the following corollary that 
will be useful later on. 

Corollary 7. Let a 6 C™" 1 be fixed. Then L(a, A) = K m X? + ™(l + o(l)) as X tends to 
infinity in the sector | argA| < it — S, and hence 

(11) Re (L(a, A)) = K m cos (^-^ arg(A)^) |A|^(1 + o(l)) 



2m 

as A — > oo in the sector | arg(A)| < n — 5. 

In particular, Re (L(a, A)) — ► +oo as A — > oo «'n any closed subsector of the sector 
|arg(A)| < a ddition, Re (L(a, A)) — > — oo as X — > oo m any closed subsector of 

the sectors < | arg(A)| < % — 5. 

Proof. This asymptotic expansion will be clear from Lemma |H1 below, or alternatively, see 
|18| Theorem 19.1] for a proof. □ 

Based on the above Corollary, Sibuya ^Sl Theorem 29.1] also proved the following asymp- 
totic expansion of the eigenvalues. 

2m 

„ / f-2n+lW \ m+2 



where 



oj = exp 



2m 



m + 2 

Notice that in this paper we consider the boundary conditions of the scaled equation (JHJ) 
where v decays in S-i USi, while Sibuya studies equation © with boundary conditions such 
that v decays in ^ U S2. The factor u m in our formula (fT2"|) is due to this scaling of the 
problem. 

Remark. Throughout this paper, we will deal with numbers like {uj u X) s for some s G R, 
and i/eC. As usual, we will use 

2-ki 



exp 



m 



and if arg(A) is specified, then 

arg {{u u X) s ) = s [arg(a/) + arg(A)] = s 



ReM 



2tt 



+ arg (A) 



m + 2 

If s 1i then the branch of X s is chosen to be the negative real axis. 



s G 



Next, we provide an improved asymptotic expansion of L. We will use this new asymptotic 
expansion of L to improve the asymptotic expansion (|12j) of the eigenvalues. 



Lemma 8. Let m > 3 and a G C m be fixed. Then there exist constants K m j(a) G C, 
< j < f + 1, such that 

L(a,X) = 



m + l 

E ? =o K mJ (a)Xi 



+ Oi\X\ 



1 

2 177 



if m is odd, 



+ 1 / \ \ 1 [ !~J 



E/=o K md (a)X2- 
as A — > 00 in the sector \ arg(A)| < % — 5. 



6m 1 n (a) i 

^HX) + |A| 



if m is even, 



Proof. The function L(a, A) is defined as an integral over < t < +00 in Theorem |H1 We 
will rotate the contour of integration using Cauchy's integral formula. In doing so, we need 
to justify that the integrand in the definition of L(a, A) is analytic in some domain in the 
complex plane. 



Let < 5 < 



m+2 



be a fixed number. Suppose that < arg(A) < n — 5. Then if 



< arg(t) < — arg(A), there exists M > such that 

6 6 5 

-7T < -- < arg(r + P m -i(t)) < arg(A) + - < tt - -, 

provided that \t\ > M . Since t m + P m _i(t) lies in a large disk centered at the origin for 
\t\ < Mo, we see that for all A with |A| large, we have that — | < arg(t m + P m _i(t) + A) < n— | 
and \t m + P m _i(t) + A| > for all t in the sector < arg(t) < ^arg(A), and hence 
^t m + P m _i(t) + A is analytic in the sector < arg(t) < — arg(A) if A lies outside a large 
disk and in the sector < arg(A) < 7r — 5. 
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Let 



m + l 



't m + P m -i(t) +\ — t~2— J2j=i bj(a)t^- j if m is odd, 



Q(t,a,X) ^ >t m + P m -i(t) + X-tf- YJ =1 bjWtf-i - if mis even. 

Then, since \Q(t, a, A)| =0 (|^|~^) as t tends to infinity in the sector < arg(£) < arg(A), 
we have by Cauchy's integral formula, upon substituting t = A^r for all A with |A| large 
enough, 

i r+°° 

Q(t, a, X) dt = X™ / Q(X™r,a, X) dr, 




n / \ — \ m-j-1 ttl „ 

P — 1 f mT) - rf - 6,-(a)^ ) if m is odd, 

A™ 



. i 



f m i(a m T ) _ r ™ _ ^ a 6-(a) T ^j J jj" 4 " 1 *" - ) if m is even. 



Similarly, (|TH|) holds for — 7r + 5 < arg(A) < 0. 

Next, we examine the following square root in Q(X™t, a, A): 



A V A(r m + 1) 



fc= 

<x> 



-Pm-l(A m T) 

^ l X(r m + 1) 



j=l Am 

where <7j(t) are functions such that #j(t) are all integrable on [0, R] for any R > 0. Moreover, 
by the definition of bj in (jHJ), we see that for 1 < j < m — 1, 

9j( T ) = / ~ r f° r some constants bj ^(o) such that > &, fc(a) = frj(a). 

£^(T™ + l) fc ~ tt ' ' 



Thus, 



T 

T 2 



t^. " ' V(^ m + i) fc ' 



L,_1 \ / 



fc=l 



777 -L 1 

O ( -^Z7 ) for all 1 < j < — ^- 



t— >oo \ T 2 



So J °° |^-(t) — bj(a)r ™ J ' dr < +00 for all 1 < j < Next, when m is even and 

= y + 1, we write 

% +i(a) \ 

#f + i(r) — x dr 

V r + A ™ / 

; P + 6™.+i(a) 

= -^m,f +i(«) ^— ln(A), 

where we take Im (ln(A)) = arg(A) G (— 7r, 7r). 
Thus, we have that 

m + l 



L(a, A) 



X^Jo -^m,j(a-)A2 + ~ + O ( I A I ~ J if m is odd, 

Tlij=Q K m j(a)X^ + ~^ - ln(A) + O ^|A|~m j if m is even, 

as A — > 00 in the sector | arg(A)| < n — S, where 

K mfi (a) = K m = J (Vl + t m - Vt^j dt>0 for all m > 3, 
(14) K m>j (a) = J ( 9j (t) - b 3 (a)t^) dt for all l<j< 



fr-+i( a ) . 

K mt rn +X {a) = j ( * ) rft when m is even. 



This completes the proof. □ 

3. Eigenvalues are zeros of an entire function 

In this section, we will prove that the eigenvalues are zeros of an entire function. 
First, we let 

G k (a) := {uj- k a 1 ,u- 2k a 2 , w"^" 1 )* a m _x) for k G Z. 

Then recall that the function f(z, a, A) in Theorem |H1 solves (jOJ) and decays to zero exponen- 
tially as z —>■ 00 in So, and blows up in S-i U Si. Next, one can check that the function 

f k (z,a,\) := f(LU- k z,G k (a),u- mk \), 

which is obtained by scaling f(z, G k (a), oj~ mk X) in the z-variable, also solves ©. It is clear 
that fo(z, a, A) = f(z, a, A). Also, fk(z, a, A) decays in Sk and blows up in Sk-i U Sk+i since 
f (z , G k (a) , uj~ mk X) decays in Sq. Since no nonconstant solution decays in two consecutive 
Stokes sectors (see LemmaEl(ii)), fk and fk+i are linearly independent and hence any solution 
of (JUJ) can be expressed as a linear combination of these two. Especially, there exist some 
coefficients C(a, A) and C(a, A) such that 

(15) f-i{z, a, A) = C(a, X)f {z, a, A) + C(a, X)fi{z, a, A). 
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We then see that 

n( u W_i,i(a, A) ~ W-i, (a, A) 

16 C a, A ) = — — — — and C a, A = - 

Wo,i(a,A) VK 0l i(a, A) 

where Wj jk = fjfa — f'jfk is the Wronskian of fa and fa. Since both fa and are solutions 

of the same linear equation (JHJ), we know that the Wronskians are constant functions of z. 

Also, since fa and fa + \ are linearly independent, W k) k+i 7^ for all k G Z. Moreover, we 

have the following lemma that is useful later on. 

Lemma 9. Suppose k, j G Z. Tnen 

(17) M4 + ij+i(a, A) = u;- 1 W fc , i (G(a),a; 2 A), 

and W ,i(a, A) = 2u^ a \ wh ere 

, s ( ¥ if m is odd, 

Moreover, 



1 



4 

— 6™ + i(a) z/m zs even. 



~ _ WLi,o(a,A) _ ^(G-Ha),^) _ 1+Ma) 
C(a ' A) " _ W ,i(a,\) ~ ~ U " ~ U ■ 

where 

if m is odd, 



18 i/a 



_i(a) i/ m even. 



Proof. See Sibuya [T%[ pages 116-118] for proof. Here, we mention that by (JHJ), we have 
6™ + i(G _1 (a)) = — 6m +1 (a) and hence y(G _1 (a)) = — z/(a). □ 



Now we can identify the eigenvalues of if as the zeros of the entire function A i— ► C(a, A). 

Theorem 10. For eacn /zxerf a G C m-1 , tae function A i— > C(a, A) entire. Moreover, A 
an eigenvalue of H if and only if C(a, A) = 0. 

Proof. Since W 0: i(a,X) ^ and since H / _i i i(a, A) is a Wronskian of two entire functions, it 
is clear from ()16|) that C(a, A) is an entire function of A for each fixed a G C m_1 . 

Next, suppose that A is an eigenvalue of H with a corresponding eigenfunction u, then 
the scaled eigenfunction v(z, A) = u(—iz, A) solves (0) and decays in S-i U Si. Hence, v is 
a (nonzero) constant multiple of fa since both decays in S\. Similarly, v is also a constant 
multiple of Thus, /_i is a constant multiple of fa, implying C(a, A) = 0. 

Conversely, if C(a, A) = 0, then fax is a constant multiple of /i, and hence fa also decays 
in S-i. Thus, /i decays in S-i U Si and is a scaled eigenfunction with the eigenvalue A. □ 

Moreover, the following is an easy consequence of (fT5Jl : For each fceZwe have 

(19) W- 1>k {a, A) = C(a, X)W 0>k (a, A) + C(a)W ltk (a, A), 

where we use C(a) for C(a, A) since it is independent of A. 



A^oo 



A— >oo 



• 2j-m-2 2j-m-2 

a) J u; 4 4 4- (1 



li 

4. Asymptotic expansions of C(a, A) 

In this section, we provide asymptotic expansions of the entire function C(a, A) as A —>■ oo 
along all possible rays to infinity in the complex plane. 

First, we provide an asymptotic expansion of the Wronskian of /o and fj in preparation 
for providing an asymptotic expansion of C(a, A). 

Lemma 11. Suppose that 1 < j < y + 1. Tnen /or eac/i a G C m_1 ; 

(20) W j(a, A) = [2iu~i + o(l)] exp [L(G J (a), u; 2j '- m - 2 A) + L(a, A)] , 
as A — ► oo in the sector 

(21) -7T + 5 < n - + 5< arg(A) < vr - 5. 
v 7 _ m + 2 ~ w ~ 

Proo/. We fix 1 < j < f + 1. Then, 

W j(a, A) = / (z, a, A)/j(z, a, A) - a, X)fj(z, a, A) 

= aT*/(0, a, A)/'(0, G^(a), ^- m - 2 A) - /'(0, a, A)/(0, (P{a)^^\) 

exp [L(G J (a),^- m - 2 A) + L(a, A)] 

[2^-2 + o{l)] exp [L(G J '(a),a; 2i - m - 2 A) + L(a, A)] , 

where we used (J§J) and (fTUj) with 

| arg(A) | < 7T - 8 and | arg(o; 2j '- m - 2 A) | < tt - 8, 

which is, ($3$. Here we also used j < — + 1. □ 

Next, we provide an asymptotic expansion of W-±,i(a, A) as A —>■ oo along the rays near the 
negative real axis. Notice from (fTBj) that V4 / _i i i(a, A) = Wo,i(a, A)C(a, A). Also, W ,i( a ; A) is 
a nonzero constant function of A. So from these one gets an asymptotic expansion of C(a, A). 

Theorem 12. For each fixed a G C m_1 and < 5 < 

(22) W_i,i(a,A) = [2i + o(l)]exp [L(G _1 (a), c^ 2 A) + L(G(a), uT m A)] , 
as A — > oo a/ono tae rays in the sector 

(23) vr - + 5 < arg(A) < tt + - o. 

m + 2 m + 2 

Moreover, there exists a constant Mi > such that W_ 1(1 (a, A) ^ for all A in the sector 
fSJ i/|A| > Mi. 

Proof. This is an easy consequence of Lemma fTTI and equation ()17|) . 

The last assertion of the theorem is a consequence of the asymptotic expansion ()22j) . □ 



The asymptotic expansion of C(a, A) in a sector near the positive real axis is obtained in 
the following theorem. 
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Theorem 13. Suppose that m > 4. Then for each fixed a G C m_1 and < 5 < ^ry 
C( a ,A) =[oj5 + (l)]exp [L(G _1 (a), cj _2 A) - L(a, A)] 

+ [cjh+Ma) + ( 1 )] exp [L(G(a), u; 2 A) - L(a, A)] , 

as A — ► oo in the sector 

41 -k An 

(24) 7T - + 5 < arg(A) < tt - — 5. 

v; m + 2 - sw - m + 2 

Proof. Suppose 2 < < y. Then from (JTTJl . (fl9|) and Lemma ITT| 

V^_i |fe (q,A) ^W^M 
u^a, aj — —— — -. — — U{a — — - — — ■ 

W 0>k {a,X) Wb,*(M) 

_ujW^ k+1 {G-\a),uo- 2 \) ~ ^- 1 W ^ 1 (G(a),uj 2 X) 



[W3 + 0(1)] 



Wo,fc(a,A) W ,fc(a,A) 

exp [L(G fc (a),^ m - 2 A) + £(£-», c^ 2 A)] 



exp [L(G fc (a), cu 2fc - m ~ 2 A) + L(o, A)] 

exp [L{G k {a),u 2k - m - 2 X) + L{G{a),u 2 X)] 



[u * + o{l))C{a, exp[L{Gk{aluJ 2k- m -2 X)+L{ajX)] 
=[ws + o(l)]exp [^((7 _1 (a), a; _2 A) - L(a, A)] 

+ + o(l)]u l+2v{a) exp [L(G(a), a; 2 A) - L(a, A)] , 

as A — > oo such that 

4(fc + 1)tt . ^ , _ 2 .. ^ 

-7r < 7r \-5 < arg(u; A) < 7T - 5, 

m + 2 

Akir 

it — + 6 < arg(A) < tt - o, 

m + 2 

A(k - 1)tt r ^ , 2u ^ 

tt i + 5 < arg(o; 2 A) < tt - S, 

m + 2 

that is, 

4&7T 47T 

tt — + 5 < arg(A) < tt — - 5, 

m+2 m+2 

provided that 2 < fc < y. So in order to complete the proof, we choose 



m— 1 



if m is odd and m > 5, 



A; = < 2 

y if m is even 



□ 



The sectors f!2Hj) and (|24jl do not cover the entire complex plane near infinity. The next 
theorem covers a sector in the upper half plane, connecting the sectors (J23J) and (|24j) in the 
upper half plane. 

Theorem 14. Suppose that a E C m_1 and < 5 < ^r^- If m>A then 

C(a, A) = + o(l)l exp [L(G _1 (a), w _2 A) - L(a, A)] 
(25) - ^!+^)+4,(a) + ^)] exp [- L (G 2 (a),u 2 - m \) - L(a, A)] , 
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as X — > oo in the sector 



(26) TV — + 5 < arg(A) < n - S. 

m + 2 



If m = 3 then 



C{a, A) =[-cj~ 2 + o(l)] exp [L(G 4 (a), cu~ 2 A) - L(a, A)] 

- [iuji +o(l)]exp [-L(G 2 (a),a;~ 1 A) - L(a, X)] , 

as A — > oo in the sector 

(27) + 5 < aig(X) < n - 5. 

5 

Moreover, if m > 6 i/ien t/iere exists a constant M 2 > st/c/i i/iai C(a, A) 7^ /or a// A in 
toe sector (J2BJ) if |A| > M 2 . 



Proof. Suppose that m > 4. Then from (|17p. (fTU|) and Lemma fTTl 

r(n \\ _ W-i,2(a,A) ^ v ^(o, A) 
C(a ' A) " Wb, 2 (a,A) " C(a) W) 

_u;^ , 3 (G'- 1 (a),^ 2 A) _ ~ u^Wo^a), a; 2 A) 



W , 2 (a,A) v y Wb, 2 (a,A) 

w[2icj-i + o(l)]exp[L(G' 2 (a),cu 2 - m A) + L{G-\a),uj- 2 X)} 
[2iuH + o(l)}exp[L{G 2 (a),io 2 - m X) + L(a, A)] 
^ _1 H / ,i(G'(a),cu 2 A) 



C(a 



[2^ + o(l)] exp [L(G 2 (a), u; 2 " m A) + L(a, A)] 
= M + exp [£(G _1 (a), uj- 2 X) - L{a, A)] 

_ a; l+2 i ,(a) 2a; M(G(a)) 

~ [2i + o(l)] exp [L{G 2 {a), u 2 ~ m X) + L{a, A)] 

as A — ► 00 such that 

—TV + 5 < 7r h (5 < arg(u;~ 2 A) < n — 5 and n h 5 < arg(A) < 7r — S, 

m + 2 m + 2 



that is, 



87T r ... r 

7r h < arg(A) < 7r — 0. 

m + 2 



Next, we use 2u(a) + /x(G(a)) = /i(a) + 4z/(a) to get ()25|). 
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Suppose m — 3. Then u b = 1. Also, M / _ 3]0 (a, A) = ^^(a, A) since /_ 3 (z,a, A) 
f 2 (z,a,X). Thus, we have that 

_ W-i, 2 (g, A) ~ Wi, 2 (a,A) 

_ g;-W_3,o(G 2 (a)^ 4 A) ~ ^T%W ^ 
W , 2 (a,A) W W 0i2 (a,X) 

u- 2 W 0>2 (G 2 (a),uj*\) ~ a^VM^aWA) 

Wo, 2 (a,A) W W ,2(a,A) 

^-^(^(a),^-^) ~ u- 1 W , 1 (G(a),u 2 \) 
W , 2 (a,X) W W , 2 (a,X) 

_ uo- 2 [2iuj-i + o(l)}exp[L(G\a),u- 2 X) + L(G 2 (a), t^A)] 
[2iw-i + o(l)]exp[L(G 2 (a),a;- 1 A) + L(a, A)] 

, u- 1 Wo, 1 (G(a),u 2 X) 

[a) [2iuj- 1 + o(l)] exp [L(G 2 (a), u^A) + L(a, A)] 
= [-aj- 2 + o(l)]exp [L(G 4 (a),a;- 2 A) -L(a,A)] 

_^l+2^(a) 2u; A t (G(a)) 

^ [2i + o(l)] exp [L(G 2 (a), u; 2 ~ m A) + L{a, A)] ' 

as A — > oo such that 

— 7r + £ < 7r — - — h<5< arg(o; _1 A) < 7r — 5 and 7r — - — \- 5 < arg(A) < n — 5, 
5 5 



that is, 



67T 

7r h 5 < arg(A) < tx — 5. 

5 



In order to show the last assertion, we suppose that C(a, A) = for some A in (I26|) with 
large |A|. Then from the asymptotic expansion ()25|). we have 

(28) exp [L{G- X (a), oo~ 2 X) + L(G 2 (a) , lu 2 '™ X)] = [ lu) -\+^)+^) + (i)]. 

By Corollary [7[ 

Re (L{G-\a),u- 2 X)) = K m cos arg(u;- 2 A)) |A|^(l + o(l)) 

= K m cos (-- + ^arg(A)) |A|* + ^(1 + o(l)), 
\ m 2m 1 

arg(^ 2 - m A)) |A|H£(i + (i)) 



Re (L(G 2 (a),cu 2 - m A)) = if m cos arg(cu 2 ~ m A)^ 



_K m cos + !!L±^ arg(A ) ) (A(-- ! --(H „(!)). 
y m 2m 
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Note that if m > 6, then < 5 < arg(A) < tt - 5 in (23). Since 

/ 2i m + 2 A /2tt m + 2 
cos I 1 arg(A) — cos I 1 arg(A) 



m 2m ) \m 2m 

/2n\ /m + 2 . . \ 

2 sin — sin arg(A) > 0, 

\ m I \ 2m I 



we see that 



Re (L(G-\a),uj' 2 X) + L{G 2 (a) , u 2 ~ m X)) -> +oo, 

as A — > oo in ()26|) . and hence the left hand side of ()28|) blows up. Thus, C(a, A) cannot have 
infinitely many zeros in (J26j) . This completes the proof. □ 

The next theorem covers a sector in the lower half plane, connecting sectors (j2SJ) and (|24jl . 

Theorem 15. Suppose that a e C m_1 and < 5 < If m > 4 then 

C(a, A) = + o(l)] exp [-L(a, cj^ 2 A) - L(G- 2 (a), cj~ 4 A)] 

+ [ U \+Ma) + ( X )] exp [L(G(a), a;- m A) - L(a, u;- m ~ 2 A)] , 

as A — > 00 in the sector 

(29) n + 5 < arg(A) < tt + — 5. 

v 7 _ v ' ~ m + 2 



If m = 3 t/ien 



C(o, A) ={-iui + o(l)] exp [-L(a, uT 5 A) - L(G~ 2 (a), cj" 4 A)] 
+ [a; 3 + o(l)] exp [L(G(a), uT 3 A) - L(a, w _5 A)] , 



as A — ► 00 in the sector 



(30) n + 5 < arg(A) < 2tt + - - 8. 

5 



Moreover, if m > 6 i/ien i/iere exists a constant M 3 > st/c/i i/iai C(a, A) 7^ for all A m 
t/ie sector (I2TH) »/[A| > M 3 . 



16 



Proof. Suppose that m > 4. Then from (fTTj) . (JT5J) and Lemma ITT1 
iy_i,- 2 (g,A) ~ Wi,- 2 (a,\) 

_WV(G- 2 (a),^- 4 A) + ~ u~ 1 W ^ 3 (G(a),u 2 X) 



~W G , 2 {G-\a),u]-*X) y J u*Wo,2(G- 2 (a),uj-*\) 

= W 0A (G- 2 (a),u- 4 \) 

[2ia;-i + o(l)] exp [L(a, cj-™- 2 A) + L(G- 2 (a), u;~ 4 A)] 

~ + o(l)] exp [L(G(a), u~ m X) + L(G- 2 (a), cu- 4 A)] 

+ o(l)] exp [L(a, u~ m ~ 2 X) + L(G- 2 (a), u^X)} 

= [2^ + o(l)] exp [L(a, a;- m " 2 A) + L(G" 2 (a), w" 4 A)] 

, r -1 , ,,1^.) exp [L(G(a),^- m A) + L(G~ 2 (a), a;' 4 A)] 
[ 1 JJ exp[L(a,w- m - 2 A) + L(G- 2 (a),cc>- 4 A)] ' 

as A — > 00 such that 

(31) 7r \- 5 < arg(c<j _4 A) < 7r — 8 and 7r — \- 5 < arg(c<j~ 4 A) < n — 5, 

m + 2 m + 2 

that is, 

(32) „ + tf< axg (A)< 7 r+ J*J--5, 

m + 2 



which is 

Suppose that m — 3. Then, 

C(a,A)=^M_c(a^ 1 '- 2(a ' A) 



W Q ,_ 2 (a,A) v 7 WV 2 (a,A) 
Wb,i(G- a (o),a;- 4 A) ~, ^^(^a) VA) 
"Wb,2(G- 2 (o),w- 4 A) W u;W , 2 (G- 2 (a),u;- 4 A) 
iyo,i(G- 2 (a),u;- 4 A) 2 ~. ^ W , 2 (G(a), cu" 3 A) 



~W , 2 (G- 2 (a),u-±\) v ; Wo, 2 (G- 2 (a),a;- 4 A) 

= W 0A (G- 2 (a),u- 4 X) 

[2iur§ + o(l)] exp [L(a, w" 5 A) + L(G~ 2 (a), u;" 4 A)] 
_ ^ 2 ~ [2jaH + o(l)] exp [L{G 3 {a),u- 4 X) + L(G(a), a;~ 3 A)] 
[2zaj-i + o(l)] exp [L(a, a;- 5 A) + L(G- 2 (a), cj- 4 A)] 

2u; M(G- 2 (a)) 

~ [2ZW" 1 + o(l)] exp [L(a, oo~ 5 X) + L(G~ 2 (a), w- 4 A)] 

, r 2 , ,l+2,(a) eX P [^(C" 2 («) ,^~ 4 A) + L(G(a) , <J~ 3 A)] 

[ 1 JJ exp[L(a,cu- 5 A) +L(G- 2 (a),a;- 4 A)] ' 

as A — > 00 such that 

7r \- S < arg(u;~ 3 A) < it — 5 and 7v — - — h 5 < arg(u;~ 4 A) < n — 5, 

5 5 
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that is, 

it + 5 < arg(A) < n H 5. 

5 

Finally, the proof of the last assertion of this theorem follows as in the proof of Theorem 

M □ 

From the asymptotic expansions in the previous four theorems, one obtains the order of 
the entire function A i— > C(a,X). The order of an entire function g is defined by 

log log M(r,g) 
hm sup , 

r^oo log r 

where M(r,g) = m&x{\g(re' l0 )\ : < 9 < 2n} for r > 0. If for some positive real numbers 
a, Ci, C2, we have exp[cir°"] < M(r,g) < exp[c2r CT ] for all large r, then the order of g is a. 

Corollary 16. The entire function X i— > C(a, A) is of order ^ + 

Proof. The sectors in ()21|) . ()23|). ()26|) . ()29|) . cover a neighborhood of infinity in the complex 

for 



2 + m 



plane. So the nonconstant entire function \C(a, A)| is bounded above by exp ci|A| 
some constant C\ > 0. Also, along the ray arg(A) = n, one can see from (fTTj) and (}2"2"|) that 
\C(a, A) | is bounded below by exp 



c 2 | A 



1 +- 

2 ' m 



for some constant C2 > 0. Hence, the order of 



C(a,.)isi + i □ 

Remark. Since the eigenvalues are the zeros of the entire function A i— > C(a, A) of order 
| H — G (0, 1), there are infinitely many discrete eigenvalues as was already mentioned in 
Theorem ^ 



5. Asymptotic expansion of the eigenvalues: Proof of Theorem El 

In this section, we prove Theorem |2] by using the asymptotic expansions of C(a, A) and 
L(a, A). 

Proof of Theorem [H Recall that by Theorem A is an eigenvalue of H if and only if 
C(a,X) = 0. 

For m > 4 and a G C m ~ 1 fixed, suppose that C(a, A) = for some A with |A| large. Then 
from the asymptotic expansion of C(a, A) in Theorem 1131 we have 

[1 + o(l)]exp [L(G{a),uj 2 \) - L{G-\a),u- 2 X)\ = -u' 2 ^ , 

and absorbing [1 + o(l)] into the exponential function then yields 

exp [L(G{a),u 2 \) - L{G-\a),u- 2 X) + o{\)] = -c^ (a) . 
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Thus, from Lemma |S] if m is odd, we infer 
In (-lu- 2 ^) 

= L{G{a),u 2 X) - L(G-\a),ij- 2 X) + o(l) 



(33) = [K m AG(a))(uj 2 X)^ -K m ,(G-\a))(iu- 2 X)^]+o(l) 

3=0 

/2tt\ , ^ Lf+1J 
= 2iK mfi sin ( — J A^+™ + ^ c mj (a)A2 +1 ^ + (l) ; 



where 

(34) c m » = tf mJ (G(a))(w 2 )*+^ - ^(G^a))^" 8 )^, 1 < j < H+l. 

Similarly, if m is even, then from LemmalHlwe have (jHSJ) with c mj (a) in for 1 < j < — , 
and 

_ x bxn +1 (a) 8ni 



c m ,™+i(o) = if m m +1 (G(a)) - K m ™ +1 {G (a)) + 



m m + 2 

where we used &m + i(G -1 (a)) = — &™ +1 (a) = &M + i(G(a)). 

Note that there exist constants M > and e > such that the function 

(35) A i-> L(G(a),uj 2 X) - L^G' 1 (a) , uj' 2 X) + o(l) 

is continuous in the region |A| > M and | arg(A) | < e. From (J33j) we then see that the 
function (|35|) maps the region |A| > M and | arg(A)| < e onto a region that contains the 
entire positive imaginary axis near infinity. 

Thus, from (J33|) we get that for every sufficiently large n G N there exists X n such that 

Lf+iJ 

" c m j(a)Xl +± ^ + o(l) = ( 2n + 1 - ) tt/. 

V m 4- 2 / 

3=1 

Thus, 

7T 



2*X m , sin {^j AI + ™ + ^ c mij (a)Al +i ^ + o(l) = ^2n + 1 - 



Lf +1 J / \ I2n + 1- ^ 

ai + -+ v CmM n , xi + ^ + o(D = [ m ; 2 2 ( 



Let 

(36) d mJ (a) 
Then 



— if 1 < j < f + 1, 



2iK m>a sin(^) 

mj( )+m £ if m is even and 7 = 2 + 1. 

2iJsT m n sinful J 2 



Lf+iJ 



(37) AI + - + £ d mj (a)Al + ^ + o(l) = _^ W+1) ,L - 

2K mi0 sm(-) 
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Introduce the decomposition A n = A 0n + Aj n , where 



2m 



i (2n+l)yr \ m+2 Ai ( „ 
A 0,n = 7777 —p2^\ and a — = I 1 ) • 



2K mj0 sin(^)y A ,„ 
Then from f)37|) we have 

1 i |2™-|-lJ 1 1— j 

A^" = (l + ^) + E d rnA^tn^ {} + jfj + ^ 

(CO r ] 

i + E( : 
fc=l ^ 



k J \Ao,n 



0(1). 



Thus, 



_ V2 + ^; Ao!n + ^ v fc J Vv, 

Lf+iJ 



E Wo,? (l + E (* + *) (£) ) + o (A„t*) 



Ai, n . 2m /- + — \ / Ai jT j 



and hence 

A ,n m + 2 ^ V & 7 VAo,n 

^f^(ge + t -)(fc)>H^) 

(38) = ^ <W a ) A o~?- 

i=i 

Thus, one concludes = A 2n + A3„, where 

2m —-L / — -L 

(39) A 2 ,n = ~ — - 2 dm,i(a)Xo t Z and A 3 , n = o (^A ™ 

Next, from along with (|H8jl we have 

A 2 ,n + A 3 ,„ + ^2 £ ( 2 t m ) (A2 ' n + A3 ' n)fc 



2m 



fc=2 

E ^(a)A Q ™ E ( 

7 = 1 \fe=l ^ 



+ ^2 E d mA a )\n E I 2 fc m ) ( A ^ + A ^) ] + ° (AoT 



J 

Lf+iJ 



(40) = ~^T2 £ <W a ) A oi- 

i=i 
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Thus, 



2m \^ 



2 m 



m + 2 ^— ' V 

fc=2 v x fcO 



S A l™ A 3,r/ 



Lf+iJ 



+ ( E (* + ^) E (J) >i»>& 1 + o fAT»-± 

3=1 \fc=l V ' ^=0 ^ / 



(41) 

and hence 



2m 



Lf+iJ 



m + 2 ^ 



/ 2 + - 

A3 '" + ^T2^i v ib , 

fc=2 v x £=0 



5^ ^J A 2,n A 3,n 



Lf+iJ 



3=1 \fc=l v " 7 £=0 



oo /l |_ l-j\ 
A: 



X] (^J A 2,n A 3,n 



"( A oX™ 



2m 



Lf+iJ 



o — 9 ! O \ / 



m + 2 



(42) 



2m 
m + 2 



3=2 

Lf+iJ 



k=2 



2 J . / ™ / 

^» A ^ ( 

7=1 \fc=l ^ 



°° /l I 1-3 

2 m ] A fc 



2.« 



So we choose 

(43) 

where 

^4,n = " 



^3,n — A4 n + A 



5,ni 



2m ^m /| + ^;\ X 2 "^m,i(°0 
«m,2(a)A ™ — m A 2 - 



m + 2 
2m 



m + 2 V 2 



m + 2 



A 0,n A 2,n 



A5,n = O ^A 0i , 

Next, we replace A3 jn in (|4^j) by Upon iterating this process we get 

Ai, n 



A n — Ao, n + Ai jn — Ao, n I 1 . - 

= Ao, n (1 + A2,n + ^3,n) 

— Ao,n (1 + A2,n + A4 jn + As )n ) 



(44) 



A 



0,d 



( ^ -± ( — 

1 + ^ e,(a)A « + o (A^ •» 

y £=1 
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as n — > +00, that is, (J5J). 

Suppose that m = 3. For this case we will use the asymptotic expansion in Theorem IT41 
that is valid in (|27)1. Similarly to what we did for the case m > 4, if C(a, A) = then from 
the asymptotic expansion in Theorem we have 

[1 + o(l)] exp [L(G\a),uj- 2 X) + L(G 2 (a), a;" 1 A)] = -iw^. 

5 3 1 

Thus, since L(a, A) = i^3 ! o(a)A6 + i^ 3jl (a)Ae + (a) A® + o(l), we have 
L{G\a), u- 2 X) + L{G 2 {a), u^X) + o(l) 

= tf 3 ,o(G 4 (a)) (^- 2 A) f + X 3 ,i(G 4 (a)) (cu- 2 A) f + K 3 , 2 (G\a)) (^ 2 A)* 

+ K 3fi (G 2 (a)) (u^X) 1 + K 3A (G 2 (a)) (u^X) 1 + K 3 , 2 (G 2 (a)) (uT^)* + o(l) 

= K 3i0 (e -1 ^ + e" l t j At + c 3j i(a)At + c 3j2 (a)A^ + o(l) 

= -2iK 3fi sin ( J As + c 3j i(a)A^ + c 3>2 (a)As + o(l). 



, 3 , 

So the continuous function A 1— > L(G 4 (a), a; _2 A)+L(Cr 2 (a), u; _1 A)+o(l) maps a neighborhood 
of the positive real axis near infinity onto a neighborhood of the negative imaginary axis near 
infinity. Hence, there exist a sequence of X n near the positive real axis such that for all large 
enough positive integers n, 

(2tv\ 5 3 1 / 15 \ 

y] A£ + c 3jl (a)A£ + c 3j2 (a)A£ + o(l) = In (-iuMj = _ 2(n + 1)tt) i. 

From this result one concludes that the asymptotic expansion (jSJ) holds for m = 3 as well 
similarly to the proof for the case m > 4. □ 

6. Proof of Theorem El 

Proof of Theorem^ First, note from (jSJ) that arg(A n ) — > as n — > +00. 
Next, we have 



2m 
m+2 



An,nH 



(2n + 3)vr 
2K mfi sin(f) 

(2n + l)yr 2tt 



2m 
m + 2 



2K m , sin(f) 2K m , sin(f) 

2m 

2 \ m+2 

An n ( 1 H 

1 2n + 1 

2m 2 / 1 

Ao ^ tl + m + 22rz+l + °l^ 



2m7r ,1---- 
(m + 2)K mi0 sin(f) 



2 m 



( 45 ) = \n + — , oN ~ — T^VV 2 V 
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Thus, 

\ _\ 2m7F \5"i i n f\ri 

n-,+oo ( m + 2) A m _ sm {^) V 
and hence, |A n+ i — A n | — ► oo and arg(A n+ i — A n ) — > as n — > +00. Since arg(A n ) — > (and 
arg(A n+ i) — > 0) as n — > +00, there exists iVGN such that |A n | < |A n +i| if n > N. □ 

Remark. Here we will show that if a £ R m_1 , then e^(a) G R for all 1 < £ < y + 1 with 
e e (a) defined in (jHJ). 

From (fTljl one can see that K m ^{G~ l {a)) = K m j(G(a)). Next, suppose that a G R m_1 . If 
m > 4 then from 



zc m » = i (if m ^G(a))(u, 2 )^ - K m ,{G- l {a)){u~ 2 )^) 

= 1 (K m ,{G{a)){uu 2 )^ -tf mJ (G(a))(w2)H^) el, 1 < j < ^ + 1. 

So by (|HE1), rf mj (a) G R for all 1 < j < f + 1, and hence by flSJ), e e (a) G R for all 
1 < £ < f + 1. 

If m = 3 then one can show ei(a) G R for £ = 1, 2, using the formulas at the end of the 
proof of Theorem |21 
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